clc;clear;
q=2;
v=36;
eps=0.0001;
hs=18;
ha=0.3;
hb=2;
%% 计算锚链长度
l=[1,1,1,1,1];%钢管+钢桶长度
for jj=1:(22.05/0.105)
    l=[l,0.105];%将锚链长度追加到长度向量中
end
for mml=1200:10:3000
    [Hk,phyK] = fH(mml,q,hs);
    if abs(Hk-hs)<0.1 && abs(phyK(5)-85/180*pi)<0.01 && abs(phyK(215)-16/180*pi)<0.01
        H=Hk;
        jd=phyK;
        ans=mml;
    end
end